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A  STUDY  OF  HEAT  TRANSFER  IN  A  CIRCULAR  CHANNEL 


Kh.  A.  Barlybayev  and  S.  V.  Bukhnan 

Flow  of  a  heat-transfer  agent  in  circular  channels  formed  by 
two  concentric  cylinders  Is  widely  used  In  modern  heat-transfer 
apparatus.  As  a  rule,  flow  of  a  fluid  In  such  chimnels  Is  turbulent. 

The  analytical  solution  of  the  problem  of  heat  transfer  In  a  laminar 

regime,  however.  Is  of  considerable  Interest.  As  In  other  cases, 
such  a  solutloiv  which  rests  on  physically  strict  equations,  allows 
the  mechanism  of  the  phenomenon  to  be  revealed  and  a  number  of  quali¬ 
tative  laws  which  are  preserved  In  turbulent  flow  to  be  established. 
Within  the  framework  of  the  problem  In  question,  let  us  first  of  all 
consider  the  effect  of  the  geometry  of  the  channel  (the  ratio  of  the 

diameters  of  the  outer  and  Inner  pipes  forming  the  circular  channel), 

which  does  not  reduce,  naturally,  merely  to  the  Introduction  of  an 
equivalent  diameter  Into  the  formula  for  calculation  of  heat  transfer. 

Statement  of  the  Problem 

Let  us  examine  the  steady-state,  hydrodynamlcally  and  thermally 
stabilized  flow  of  an  Incompressible  fluid  with  constant  physical 
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parameters  along  a  circular  channel  formed  by  two  coaxial  cylindrical 
pipes  (the  diagram  and  position  of  the  X  and  r  axes  are  shown  In 
Fig.  1).  The  radii  of  the  Inner  and  outer  pipes  will  be  denoted  by 
and  Rgf  respectively. 


Let  the  circular  space  between  the  Inner  pipe  and  the  cylindrical 
surface  of  radius  (Fig.  l)  on  which  the  maximum  temperatures  lie 
be  section  1,  and  let  the  space  between  the  outer  pipe  and  this  sane 
cylindrical  surface  be  section  2.  Let  us  Isolate  In  these  sections 
circular  layers  with  radius  r  and  thickness  dr.  The  difference 
between  the  amounts  of  heat  Introduced  and  removed  by  a  fluid  In  a 
section  of  length  dx  In  the  region  from  r^^^  to  r  let  us  call  Q.  Since 
on  the  cylindrical  surface  of  radius  r^^^  the  derivative  of  the  tem¬ 
perature  Is  equal  to  zero  (•^— )^  =*  0,  then  a  quantity  of  heat  Q  will 
be  transmitted  through  the  layer  of  thickness  dr.  In  this  case  the 
following  equation  holds: 

j*  ^  l/(r)  rdr  -  2k)v  5,  .  ^ 

where  p,  Cp  and  X  are  the  density,  specific  heat  and  thermal  conduc¬ 
tivity  of  the  fluid,  and  U(r)  Is  the  flow  rate  of  the  fluid. 

With  steady-state  motion,  from  the  condition  of  constancy  of 
heat  flows  through  the  walls  of  a  channel  it  follows  that  the  tem- 
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perftture  of  the  walls  of  the  channel  and  the  temperature  of  the  fluid 
vary  along  the  channels  according  to  a  linear  law.  Ihereforef  the 
change  In  temperature  of  the  fluid  In  a  channel  of  length  dx  can  be 
determined  from  the  heat  balance  for  the  entire  stream 


Hence 


« IR\  -  R\)  U^,  ^  &  («,  +  q,  Rtl 

^  _  Xti Si  + 


(2) 


where  tJ  Is  the  average  velocity  of  the  fluid  In  the  c.iannel,  and 
q2  the  heat  flows  through  the  Inner  and  outer  walls  of  the  channel.. 

Substituting  the  value  of  Into  (1),  we  obtain  an  equation 
which  allows  us  to  find  the  temperature  field 


I 


2  (^1  "4“  ^  Aff) 


U  ir)’rdT  !=■  rl  ^ , 


(5) 


The  radius  of  the  maximum  of  the  temperature  profile  r_  can  be 

m 

found  from  the  condition  that  the  heat  Introduced  by  the  fluid  Into 
section  1  or  2  will  be  transmitted  to  the  Inner  and  outer  pipes 
forming  the  channel,  respectively.  This  condition  for  section  1  has 
the  form: 

f  2-:.C,(/(r)^r</r-2i:«,»,. 

(M 

Substituting  the  value  of  ^  from  (2)  Into  this  expression,  we 
obtain 


^^^-^rU{r)rdr~R,q,: 


(5) 


Expression  (5)  is  the  equation  necessary  for  determining  r^j^Cln  the 
limit  of  the  Integral).  In  order  to  solve  It,  It  Is  necessary  to  know 
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the  dependence  of  velocity  upon  the  redlue 


Fluid  Plow  with  a  Constant  Velocity 

Let  us  consider  the  solution  of  the  problem  for  the  case  of  a 
uniform  velocity  distribution  along  the  cross  section  of  the  channel 
(U(r)  ■  IT  •  const).  Under  this  assumption,  the  solutions  of  Eq.  (4) 
separately  for  sections  1  and  2  have  the  form: 


*1  a  t  -  ti^  -  7*+  -P). 


(6) 


*2  s  t 


‘'2W 


(7) 


where  and  tg^  are  the  temperatures  of  the  Inner  and  outer  pipes, 
respectively.  Here  we  use  the  boundary  conditions  t  »  t^^^  at 
r  ■  and  t  «  tg^^  at  r  «  Rg. 

The  temperatures  of  the  walls  of  the  Inner  and  outer  pipes  are 
interlinked  by  the  i*elationshlp 


% 


^Iw 


4 


(8) 


If  In  Eq.  (6)  the  temperature  of  the  inner  pipe  (tj^^)  is  expressed 
in  teirnis  of  the  temperature  of  the  outer  pipe  (tg^),  according  to 
Eq.  (8),  then  for  the  temperature  at  any  point  In  the  channel  we  have 
the  equation 

t 

^  ~  ^2W  _  2r»,lp  + 

^ra  “  ^2w  2  1"-^- 

where  t  Is  the  maxlmvun  temperature,  corresponding  to  the  radius  r  . 
m  m 

The  value  of  r^^^  entering  Into  this  equation  Is  foxuid  from  (5) 
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^  ™  5'*  * 


i  +  ♦«  *11 


(10) 


Under  the  condition  of  constancy  of  heat  flows  through  the  Inner  and 
outer  walls  of  the  channel >  the  value  of  r^^^  Is  defined  as  the  geo¬ 
metric  mean  of  the  radii  of  the  tubes  forming  the  channel,  and  has 
the  form: 

VlTT?r 


Figure  2  shows  the  dlmenslonless-temperatvire  profiles  constructed 
according  to  Eq.  (9)  for  a  channel  with  Rg^  ■  2  at  various  ratios  of 
heat  flows  through  the  Inner  and  outer  pipes,  and  also  temperature 
profiles  for  various  circular  channels  (Rg^  ■  2,  10,  lOO)  at  the  same 
heat  flows  through  the  walls  of  the  channel  (qg^^  ■  l).  From  the 
figure  It  Is  apparent  that  at  qgj^  «  1  the  temperature  profile  In  a 
circular  channel  Is  unsymmetrlcal:  the  temperature  of  the  Inner  wall 
Is  higher  than  that  of  the  outer  wall.  Here  the  temperature  maximum 
Is  located  near  the  Inner  pipe.  Figure  5  shows  graphs  of  position 
of  the  temperature  maximum  of  the  fluid  as  a  function  of  the  ratio  of 
the  radii  of  the  pipes  forming  the  circular  channel  for  various  ratios 
of  heat  flows  through  the  walls.  Figure  4  gives  a  graph  of  the 
dimensionless  difference  between  wall  temperatures  as  a  function  of 
R21  for  qgj^  -  1  and  qg^^  -  2.  As  Is  apparent  from  these  figures,  the 
asymmetry  of  the  temperature  profile  and  the  difference  between  the 
temperatures  of  the  walls  Increase  with  an  Increase  In  the  ratio  1^1. 
The  location  of  the  maximum  temperatures  as  functions  of  the  ratio  of 
heat  flows  through  the  walls  are  shown  In  Figs.  5  and  6  (the  latter 

*  R  and  £  with  two  subscripts  denotes  the  ratio  of  two  dimensional 
radii  and  heat  flows;  the  first  subscript  pertains  to  the  numerator, 
the  second  to  the  denominator. 
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for  R2||^  ■■2).  It  l8  apparent  from  these  figures  that  the  location  of 
the  temperatiire  maximum  and  also  the  dimensionless  difference  between 
the  wall  temperatures  of  the  fluid  vary  with  The  temperature 

maximum  approaches  the  Inner  wall  of  the  channel  as  the  ratio 
Increases  and  approaches  the  outer  wall  as  It  decreases.  In  the 
limiting  case,  at  q2^  oo  (the  Inner  wall  Is  heat-insulated)  and  at 
qg^  *  0  (the  outer  wall  Is  heat-insulated),  the  temperatxire  maxlsaim 
will  be  on  the  Inner  or  outer  wall,  respectively.  For  Rg^  2,  the 
temperatures  of  the  walls  of  the  channel  at  qgj^  ■  O.785  (Fig.  6)  will 
be  equal  to  one  another.  When  the  ratio  qg^^  increases,  the  tempera¬ 
ture  of  the  Inner  pipe  will  become  higher  than  that  of  the  outer  one, 
and  their  difference  will  Increase  with  an  Increase  In  q2i*  With  a 
decrease  In  qg^,  the  sign  of  the  difference  between  the  wall  tempera¬ 
tures  will  change  and  the  temperature  of  the  Inner  pipe  will  become 
lower  than  that  of  the  outer  one.  This  difference  Increases  with  a 
decrease  In  the  ratio  qg^^. 


'^”^2w  2.  Temperature  dlstrl- 

butlon  In  circular  channel 
with  laminar  flow  (dotted 
line)  and  with  flow  with  a 
constant  cross-sectional 
velocity  (smooth  line). 

1)  Rgi  -  2,  qgi  -  O.785J  2) 

“  2,  qgj^  "  Ij  3)  “  2, 

^21  “  ”21  "  ^21  " 

6)  R2j_  “  ^21.  “ 

Rgj^  »  100,  qgj^  -  1;  9)  velocity 
distribution  at  Rgj^  -  lOj  ' 

10)  temperatvire  distribution 
at  "  10,  ^21  * 
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Fig.  4.  Dependence  of 
difference  In  vftH  tem¬ 
peratures  upon  ^  . 


2)  Rgi  -  10. 


Fig.  6.  Dependence  of 
difference  between  wall 
temperatures  upon  ^  for 

^1 

a  channel  with  »  2. 
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Calculation  of  Heat  Transfer 


Let  U8  examine  the  calculated  heat-transfer  dependences  for  a 

uniform  velocity  distribution  (u(r)  •  V  >  const).  As  Is  known,  the 

qd 

dimensionless  heat-transfer  coefficient  (Nusselt  number)  Is  Nu  •  » 

where  «  t^^  -  t^  Is  the  average  temperature  head,  and  d^^ 

the  equivalent  diameter  of  the  channel.  Together  with  the  calculation 
of  heat  transfer  for  the  entire  channel  as  a  whole,  let  us  attempt  to 
determine  heat  transfer  separately  for  each  of  the  boundary  walls.  In 
order  to  explain  more  clearly  the  effect  of  the  curvature  of  the  walls. 
For  this  let  us  Introduce  the  average  temperature  heads,  the  heat- 
transfer  coefficients,  the  equivalent  diameters  and,  finally,  the 
Nusselt  number  not  only  for  the  entire  circular  channel, 
but  also  for  sections  1  and  2  separately.  It  is  understood  that  this 
division,  as  before.  Is  possible  because  sections  1  amd  2  aT«  separ¬ 
ated  by  an  adiabatic  surface  with  radius  r^j^.  The  average  temperature 
heads  with  respect  to  the  flow  rate  for  sections  1  and  2  are  deter¬ 
mined  as: 


rm 

f  2k  (7  *1  r^r 

_ 

*1  av  “  ’ 

r"* 

J  i*V%,r4r 
*2  av  “  _£a _ _ 


(11) 

(12) 


Substituting  here  the  values  of  and  Ag  from  (6)  and  (7),  we  obtain 

X 

A.  ■  2  (Ti  A?,  ■»<,»,)  *  r*  (l3) 

1  X (f«,  - «,«)  («*,  -  fr,)  {  2  $)•  ' 


'2  av 


2(f(ff,  +  ft  lt*m  1.  1  8  _4  1 

(T  - 2“  +  -8''^-+  •  j-  (lA) 
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Recalling  that  -  l  (^)  - 1  (7?),-,.  and  d^^  -  2  (vR^), 

dg  •■  2  (R2''^in^*  the  Nueselt  numbers  for  sections  1  and  2 

lav  (15) 


““2  “ 

2av 


a.(4»^«  ~  -4r*m  lft+ir*m  +  ^ 


(16) 


The  results  of  the  niomerlcal  calculation  for  sections  1  and  2  at 
equal  heat  flows  through  the  inner  and  outer  walls  for  a  few  values 
of  Rg^  are  given  in  Table  1. 


TABIS  1 


q  1  av 

—  A 

q  2  av 

Nu^ 

Nug 

*2  av 
*1  av 

NUg 

“2 

“1 

2 

0.126 

0.187 

6.51 

6.29 

1.48 

0.96 

0.67 

10 

0.61 

2.09 

7.04 

6.50 

3.42 

0.92 

0.29 

100 

1.59 

24.52 

11.31 

7.34 

15.42 

0.65 

O.O65 

Here  and  Og  are  the  heat-transfer  coefficients  for  the  Inner  aoid 
outer  walls  of  the  circular  channel. 

Now  let  us  calculate  heat  transfer  for  the  entire  circular 
channel.  The  avereige  temperature  head  for  the  entire  channel  section 
let  us  represent  In  the  form 


Using  (13)*  and  {17)»  let  us  determine  the  Nusselt  number  for 
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the  entire  channel 


-  a«.)(a.  -  a.r 

4  r*m  »■  («•» -H\H **i— IFT  * 


(18) 


where  q, 


av 


V'a 


iB  the  average  heat  flow  through  both  walls , 


and  ■  2(R2-Rj^)  the  equivalent  diameter  of  the  circular  channel. 

Figure  7  shows  the  Nusselt 
number  as  a  function  of  the  ratio 
of  the  radii  of  the  outer  and 


Inner  pipes.  As  Is  apparent 
from  the  figure.  In  the  limiting 
case  of  transition  to  a  plane 
channel  (Rgj^  “  1)»  the  Nusselt 
number  becomes  equal  to  about 
12.  In  the  other  limiting  case  — 


•  U  U  4 

I 

Fig.-7.  The  Nusselt  number 
versus  the  ratio 


transition  to  a  circular  pipe  (f^i  ®')“  the  value  of  Nu^^^  approaches 
8.  Both  these  limiting  cases,  naturally,  coincide  with  the  results 
of  the  corresponding  solutions  given  in  the  literature  [1]. 


The  Effect  of  the  Velocity  Distribution 


In  order  to  evaluate  the  effect  of  the  velocity  distribution 
along  the  cross  section  of  the  channel,  which  was  not  taken  into 
account  earlier,  let  us  examine  the  solution  for  laminar  flow  in  a 
circular  channel.  Let  us  limit  ourselves  to  the  particular  case  of 
equality  of  heat  flows  through  the  walls  of  the  channel  (^21  -  ^^)- 
The  velocity  distribution  of  laminar  flow  in  a  circular  channel  under 
isothermal  conditions,  as  is  known,  satisfies  the  equation  [2] 


(19) 
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The  value  of  the  radius  correspondlns  to  the  maxlmuin  velocity  cal¬ 
culated  by  this  equation  has  the  form 

=  (20) 

^e  velocity  profiles  constmicted  according  to  Formula  (I9)  are  unsym- 
metrlcal .  The  velocity  maximum  does  not  coincide  with  the  middle  of 
the  circular  channel  and,  as  the  temperatiire  maximum  In  the  previous 
calculation  at  <i2±  “  situated  near  the  Inner  pipe. 

The  average  flow  rate  of  the  fluid  Is  determined  from  the 
equation 

fi* 

^  2»U{r)rdr 

n  _  'Ll _ Ll 

^  S’)' 


(21) 


Substituting  the  value  of  U(r)  from  (19)  and  of  U  from  (21)  Into  (5), 
we  obtain  an  equation  for  determining  the  radius  of  the  maximum  of  the 
temperature  profile 


1 


In  Rn  ^ 

=  Rii  +  Rit  • 


/?  a  - 


'"ir; 

^  n 


1*^1 


(22) 


The  graphical  solution  of  this  equation  for  r^j^  shows  that  the 
locations  of  the  temperatiire  maximum  for  laminar  flow  of  a  fluid  and 
for  flow  with  a  constant  cross-sectional  velocity  coincided,  for  all 
practical  purposes.  As  an  example.  Table  2  gives  values  of  the  radii 
of  the  maximum  of  the  temperature  profile  for  various  values  of  the 
ratio  at  qgj^  »  1.  Values  of  the  maximum  velocity  profile  for 
laminar  flow  are  given  here  for  congiarison. 

It  should  be  noted  that  the  radius  of  the  maximum  of  the  tempera¬ 
ture  profile  r^ij  for  any  value  of  Rgj^  Is  less  than  the  radius  of  the 
maxlmxun  of  the  velocity  profile  as  determined  by  Eq.  (20)  (in 
Fig.  5  the  dotted  curve  shows  R^/Rg  as  a  function  of  Rgi)* 
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12 


For  clarity*  the  velocity  and  temperature  distribution  for  lam¬ 
inar  flow  In  a  circular  channel  with  ■  iO  and  ■  1  la  given  In 
the  corner  of  Pig.  2. 

The  solution  of  Eq.  (3)  talcing  (I9)  and  (20)  Into  account  gives 
the  cross-sectional  temperature  distribution  for  a  circular  channel 
In  the  form: 


where 


-« +  B  (n  +  ]  i«  4*-  + 


SBs  +  r* 
16 


-  -f 


16 


R«i-  A?, 
■  Id  f(t, 


(23) 


Curves  of  the  temperature  distribution  constructed  according  to 
Eq.  (23)  for  various  values  of  are  given  In  Pig.  2  by  dotted 
lines.  As  Is  apparent  from  this  flgvire,  the  curves  of  the  tempera- 
tiire  distribution  for  laminar  flow  and  for  flow  with  constant  cross- 
sectional  velocity  differ  from  one  another  very  little  with  regard  to 
the  location  of  the  temperature  maximum.  This  allows  it  to  be  assumed 
that  the  location  of  the  temperature  maximum  for  a  turbulent  velocity 
profile,  without  taking  Into  account  the  effect  of  turbulent  thermal 
conductivity,  remains  practically  the  same.  This  result  pertains  to 
the  practically  Important  case  of  fluids  with  a  very  low  Prandtl 
number  (liquid  metals). 

This  conclusion  cannot  be  extended  to  the  value  of  the  tempera¬ 
ture  gradient  near  the  walls  of  a  channel  and,  therefore,  to  Iseat  . 
transfer.  As  Is  apparent  from  Pig.  2,  near  both  walls  the  temperature 
gradient  In  the  case  of  U(r)  =»  TT  *  const  Is  markedly  higher  than  when 
the  laminar  velocity  profile  Is  taken  Into  account.  This  effect  of 
the  velocity  profile  upon  heat  transfer  Is  In  qualitative  agreement 
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with  that  known  for  heat  transfer  In  a  circular  pipe  and  In  a  plane 
channel.  Let  us  reaeaber  that  In  the  first  case  ■  oo,  circular 
pipe)  the  ratio  of  the  Nusaelt  numbers  for  a  constant  velocity  and 
parabolic  distribution  equals  1.85»  while  In  the  second  case 

"  1»  plane  channel)  It  equals  «  1.45.  These  data  can  be  used 
for  an  approximate  evaluation  of  the  effect  of  the  velocity  profile 
upon  heat  transfer  In  a  circular  channel.  When  necessary,  an  accurate 
calculation  can  be  made  similarly  to  the  prevloue  one,  but  It  Is  not 
given  he]?e  owing  to  Its  awkwardness.  Let  us  note  only,  for  example, 
that  at  ^  *121  **  ^  ratio  of  Nusselt  numbers  cal¬ 

culated  with  and  without  taking  the  velocity  distribution  Into  account 
are  1.55  and  1.65  respectively. 

Conclusions 

As  can  be  seen  from  the  calculation  made,  the  geometry  of  the 
channel  has  the  main  effect  on  the  temperature  field  and  heat  transfer. 
Just  as  on  the  flow  Itself.  The  asyraraeti^  of  the  temperature  profile, 
particularly  the  location  of  the  temperature  maximum,  la  a  function 
of  the  ratio  of  the  radii.  The  same  sort  of  dependence  Is  character¬ 
istic  of  the  Nusselt  nxjunber  In  stabilized  heat  transfer,  which  In 
general  Is  determined  by  two  pearameters:  the  ratio  of  the  radii  and 
the  ratio  of  the  heat  flows  through  walls  of  the  channel. 

The  determination  of  heat  transfer  separately  for  the  inner  and 
outer  walls  of  a  circular  channel  shows  that  for  the  Inner  wall  It 
Is  more  Intensive,  especially  at  high  ratios  of  the  radii  of  the 
outer  and  Inner  pipes  of  the  circular  channel. 

The  calculations  (for  the  case  of  constancy  of  heat  flows 
through  the  walls  of  the  channel)  show  that  the  location  of  the 
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temperature  maximum  Is  slightly  dependent  upon  the  velocity  profile. 
Regarding  the  effect  of  the  latter  on  heat  transfer.  It  apparently 
has  the  saune  character  as  for  heat  transfer  In  the  two  limiting  cases 
(in  a  circular  pipe  and  In  a  plane  channel). 
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MAJOR  AIR  OOMMMDB 


nr.  OoplM 


AF8C 

SCFID  1 

ASIIA  2$ 

TSBTL  $ 

TDBDP  $ 

AZDC  (AST)  1. 

BSD  (BSP)  1 

IPPTC  (PTT)  1 
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